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Abstract 

We study the discrete spectrum of the Hamiltonian H = —A + V{r) for the 
Coulomb plus power-law potential V(r) = —l/r + (3 sgn(q)r g , where j3 > 0, q > 
—2 and q ^ . We show by envelope theory that the discrete eigenvalues E n £ of 
H may be approximated by the semiclassical expression E n i(q) w min r> o{l/r 2 — 
+ sgri(q)f3{yr) q }. Values of fx and v are prescribed which yield upper and 
lower bounds. Accurate upper bounds are also obtained by use of a trial function of 
the form, ip(r) = r i + 1 e -( xr ) q , "W e give detailed results for V(r) = —l/r + /3r q , q = 
0.5, 1, 2 for n = 1, £ = 0,1,2, along with comparison eigenvalues found by direct 
numerical methods. 
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1. Introduction 

In this paper we derive upper and lower bound formulas for the spectrum of a 
single particle in three dimensions that obeys non-relativistic quantum mechanics 
and has Hamiltonian 

H = -ujA - A/r + Bsgn(q)r q , u, A, B > 0, and q ^ 0, q > -2. (1.1) 

The Coulomb plus power-law potential is of interest in particle physics where it 
serves as a non-relativistic model for principle part of the quark-quark interac- 
tion. This class of potentials has been well studied and much work has been done 
to approximate the eigenvalues, with or without the Coulomb term necessitated 
by QCD [1-15]. Our goal in this paper is to provide simple formulas for upper 
and lower energy bounds for this class of potentials. Firstly, we use the 'envelope 
method' [16,17] to obtain upper and lower bound formulas for all the discrete eigen- 
values. We also use a Gaussian trial function and the 'sum approximation' [18,19] 
to improve the bounds for the bottom of each angular-momentum subspace. The 
energy bounds so far discussed may all be expressed in terms of the following semi- 
classical energy formula: 

L\-- + Bsgn(q)(urA (1.2) 



Rj mm 

r>0 



for suitable choices of the parameters \x > and v > 0. We also apply a variational 
method used earlier [20] which is based on the exact Coulomb wave function and 
yields accurate upper bounds for the bottom of each angular momentum subspace. 
We compare all these results with 'exact' eigenvalues computed by direct numerical 
integration. 

For the class of potentials studied some exactly solvable cases exist for suitable 
values of the couplings u, A, B, and the power q. For example, for the well- 
known hydrogenic atom and the harmonic oscillator potentials we have explicitly 
for n = 1, 2, 3, ... 

« = =* E -'--j^w <L3) 

and 

q = 2 =>- E n £ = (ujB)^(An + 2^ — 1). (1.4) 

For £ = 0, exact solutions are also available for the linear potential q = 1. We 
can simplify the coupling problem in general by the use of scaling arguments. If, for 
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each fixed g, we denote the eigenvalues of H = —ujA — A/r + Br q by £(u>, A, B) , 
and consider a scale change of the form s = r jo , and choose the scale a = u/A, 
then it is straightforward to show that, 

S(u,A,B)= (^pj£(l,l,p), /3=(?^(^y + \ (1.5) 
Hence, the full problem is now reduced to the simpler one-parameter problem 

H = -A - 1/r + (3 sga(q)r q , E = E(p) = E(l, 1, P), P > 0. (1.6) 



2. Energy bounds by the envelope method and the sum approximation 

The comparison theorem tells us that an ordering between potentials implies 
an ordering between the corresponding eigenvalues. The 'envelope method' [16, 17] 
is based on this theorem and establishes upper and lower bound formulas for a 
wide class of attractive spherically-symmetric potentials. We need a solvable model 
—A + h(r) which provides an 'envelope basis' for the study of the problems of the 
form —A + g(h(r)), where the transformation function g is monotone increasing 
and of definite convexity: when g is convex, we obtain lower bounds; when g is 
concave, the theory yields upper bounds. The natural basis in this contex is a single 
power-law potential. The spectrum of a Hamiltonian of the form 

H = —A + sgn(g)r 9 , where q > — 2 and g / (2-1) 

may be represented exactly by the following semiclassical expression [12, 17], 

E ne = mini ^ + sga(q){P ne (q)r) q \ (2.2a) 

= sg „(,)( 1 + |)(^)*. (2.2,) 

The function P = P n e(q) is known as the P-representation, for the Schrodinger 
spectra generated by the power-law potentials. It is convenient to use the P function 
to study and analyse the spectra of these problems mainly because it is known [12] 
that P n e(q) is monotone in q and it is also smoother than E n £ as a function of 
g; the case g = corresponds exactly to the log potential. From (1.3) and (1.4) 
we find, in particular, that: 

P ne (-l)=n + £ (2.3) 
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and 

P ne (2) = 2n + £+l/2. (2.4) 

In Table 1 we exhibit some numerical values for P n i{\) and P n g(l). We have found 
the exact eigenvalues for the linear potential in terms of the zeros of the Airy func- 
tion, but those for q = | have to be computed numerically: this use of some isolated 
numerical input is justified since, for each {n, £} pair, the resulting approximation 
formulas include all the potential parameters but depend only on a single 'numeri- 
cal input'. Envelope theory [13, 18] shows that the eigenvalues of the Coulomb plus 
power-law potential may be approximated by the following semiclassical expression, 

f 1 1 ] 

£ pa min < — h (3 sgn( q)(ur) q > , where u, v > 0. (2.5) 

Since V(r) = g(h(r)) is at once a convex function of h(r) = — 1/r and a concave 
function of h(r) = sgn(g)r 9 , the spectral representation P n i(q) allows us to specify 
upper and lower bounds formulas as follows. If fj, = v = P n e(— 1) , then £ is a 
lower bound for E n g , and if = v = P n i(q), then £ is an upper bound. We 
may improve the lower bound for the bottom of each angular momentum subspace 
by using the sum approximation [18,19], which is equivalent to the choice \i = 
Pu(— 1) = (£ + 1) and v = Pu(q). For the bottom of the spectrum we can also 
improve the upper bound by using a Gaussian trial function and minimizing over 
scale: this is equivalent [13] to using the parameter values 

=**=(§)* 

We note that the same parameters and v which guarantee that (2.5) yields vari- 
ous energy bounds may also be used in the 'full' semiclassical formula (1.2), includ- 
ing all the original Hamiltonian parameters {cu, A, B}. In Section 3 we apply (2.5) 
to the explicit cases V(r) = — 1/r + (3 r q for i = 0, 1, 2, where q = 1, 2, and 0.5 . 



2r((3 + <z)/2) 



(2.6) 



3. Variational method 

The second approach in this paper is to use a trial function explored in pre- 
vious work [20] to obtain accurate upper bounds for the bottom of each angular 
momentom subspace. We start with Schrodinger's equation 

#V(r) = (-A --+(3 8gn(q)r*) ^(r) = E n£ ([3)ij(r), q ^ 0, q > -2. (3.1) 
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This problem is solvable if (3 = 0, and the corresponding wave function ifj(r) is 
given by 

^(r) = r i+1 e- xr L 2 J +1 (2xr). (3.2) 

In order to obtain an upper bound for the bottom of each angular momentum 
subspace E\i for fixed power q we choose ip(r) to be of the following form 

^(r) = r e+1 e-^ d (3.5) 

and define S by S(/3,x,d) = ^j^jy-, where x and d are variational parameters. 
Now, we minimize £ with respect to x and d. The necessary conditions for a 
critical point are ff = and ff = 0. Consequently, using (3.1) and (3.5), we 
obtain the following upper bound formula for the eigenvalues En 

£u(/3 7 d, x) = a\x 2 — a,2X + asx~ q , (3.6) 

where a± , ci2 and 03 are as given below 

2-2 

dl = 2 d 



i 2^(2£ + l)(2£ + rf+l)r(^i) 



fl 2^^ 
a 2 = z d ■ 



a 3 = sgn(g)/32^r - ;/ 



By using (3.6) we derive the following equation for x 

x ^ _ ^L^+i -^3=0. (3.7) 
2ai 2ai v ' 

After solving Eq.(3.7) to obtain x, from the numerical solution of ^| = we find 
<i for n = 1 and £ = 0, and then we use the same d value for all i. 



4. Results and Conclusion 

We have found general semiclassical energy formulas (1.2) and (2.5) for the 
eigenvalues generated by the Coulomb plus power-law potentials. Specific values for 
the parameters fx and v are given which guarantee that the formulas yield bounds 
for all the discrete energies. By using a more finely tuned wavefunction, we have 
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also derived an improved upper bound (3.6) valid for the bottom of each angular 
momentum subspace. We may rewrite (2.5) in the form of a pair of parametric 
equations for the curve {/?, E(f3)} . For fixed q > — 1 we obtain: 



By envelope theory, we know that these parametric equations yield a lower bound 
if n = v = P n g(—1) = (n + £), and an upper bound when fj, = v = P n g(q). For 
the bottom of each angular momentum subspace the prescription \i = Pu(— 1) = 
(£ + 1), v = Pu(q) yields an improved lower bound. An improved upper bound 
for the bottom of the spectrum is given by using the 'Gaussian' P-numbers (2.6). 
In Figures 1, 2, and 3, we plot the function E(f3) for n = 1 ,£ = 0, 1, 2 for the 
Coulomb plus harmonic oscillator (q = 2) , Coulomb plus linear (q = 1), and 
Coulomb plus r 5 potentials, along with the corresponding accurate variational 
bounds using (3.6) (dashed line), and some comparison numerical values represented 
as stars. The advantage of the semiclassical formulas is that that they describe in 
approximate analytical form how the eigenvalues depend on all the parameters of 
the problem. 
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Table 1 The 'input' values P n i{^) and P n t(l) to be used in the general formula 
(2.5) for the energies corresponding to the potential V(r) = — 1/r + f3 sgn(q)r q . 
These P- values yield upper bounds when q < |, or q < 1, respectively. 



n 


I 




Pnt(l) 


1 





1.30266 


1.37608 


2 





2.97387 


3.18131 


3 





4.65440 


4.99255 


4 





6.33742 


6.80514 


5 





8.02149 


8.61823 


1 


1 


2.29747 


2.37192 


2 


1 


3.93966 


4.15501 


3 


1 


5.60154 


5.95300 


4 


1 


7.27194 


7.75701 


5 


1 


8.94679 


9.56408 


1 


2 


3.29535 


3.37018 


2 


2 


4.92261 


5.14135 


3 


2 


6.57089 


6.92911 


4 


2 


8.23022 


8.72515 


5 


2 


9.89619 


10.52596 


1 


3 


4.29424 


4.36923 


2 


3 


5.91240 


6.13298 


3 


3 


7.55077 


7.91304 


4 


3 


9.20118 


9.70236 


5 


3 


10.85929 


11.49748 


1 


4 


5.29352 


5.36863 


2 


4 


6.90560 


7.12732 


3 


4 


8.53658 


8.90148 


4 


4 


10.17964 


10.68521 


5 


4 


11.83110 


12.47532 
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Figure 1. 

The eigenvalues E(f3) of the Hamiltonian H = —A — 1/r + f3r 2 for dimension 
N = 3, n = 1, and I = 0, 1, 2 . The continuous curves show the upper bound EU 
given by the envelope formula (2.5) with v = n = Pn(2), for i = 1,2 and the 
lower bound ELS by the sum approximation given by the same formula but with 
v = Pu(2) and n = Pu(— 1). The upper bound for £ = is calculated using 
v = -P^(2) and [i = P^(—l) in formula (2.5). The dashed curve EC represents 
the upper bound by formula (3.6). The stars EX represent accurate numerical data. 
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Figure 2. 

The eigenvalues E((3) of the Hamiltonian H = —A — 1/r + (3r for dimension 
N = 3, n = 1, and I = 0, 1, 2 . The continuous curves show the upper bound EU 
given by the envelope formula (2.5) with v = fx = Pn(2), and the lower bound 
ELS by the sum approximation given by the same formula but with v = Pu(l) 
and ii = P\i{—\). The upper bound for £ = is calculated using v = -P^(l) and 
li = Pie(— 1) in formula (2.5). The dashed curve EC represents the upper bound 
by formula (3.6). The stars EX represent accurate numerical data. 
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Figure 3. 

The eigenvalues E((3) of the Hamiltonian H = —A — 1/r + /3r ' 5 for dimension 
N = 3, n=l, and I = 0, 1, 2 . The continuous curves show the upper bound EU 
given by the envelope formula (2.5) with v = fx = P^(0.5), and the lower bound 
ELS by the sum approximation given by the same formula but with v = Pu(0.5) 
and ii = Pn(—1). The upper bound for £ = is calculated using v = P^(0.5) 
and ii = Pit(— 1) in formula (2.5). The dashed curve EC represents the upper 
bound by formula (3.6). The stars EX represent accurate numerical data. 
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